SEMI-CLASSICAL SCHRODINGER EQUATIONS WITH 
HARMONIC POTENTIAL AND NONLINEAR PERTURBATION 

REMI CARLES 

Abstract. Solutions of semi-classical Schrodinger equation with isotropic har- 
monic potential focus periodically in time. We study the perturbation of this 
equation by a nonlinear term. If the scaling of this perturbation is critical, 
each focus crossing is described by a nonlinear scattering operator, which is 
therefore iterated as many times as the solution passes through a focus. The 
study of this nonlinear problem is made possible by the introduction of two 
operators well adapted to Schrodinger equations with harmonic potential, and 
by suitable Strichartz inequalities. 

Resume. Les solutions de I'equation de Schrodinger semi-classique avec poten- 
tiel harmonique isotrope focalisent periodiquement en temps. Nous etudions 
la perturbation de cette equation par un terme non-lineaire. Pour une echelle 
critique de cette perturbation, chaque traversee dc foyer est dccritc par un 
operateur de diffusion non-lineaire, qui est par consequent itere autant de fois 
que la solution traverse une caustique. Cette etude est permise par I'usage de 
deux operateurs qui s'averent bien adaptes a I'equation de Schrodinger avec 
potentiel harmonique, et par des estimations de Strichartz adequates. 



1. Introduction 
Consider the initial value problem, 

where e g]0, 1] is a parameter going to zero and / is a smooth function, say / S 
iS(M). The potential is the isotropic harmonic potential, 

(1.2) Vix)^\^^lia:l + --- + xl). 

The case of anisotropic harmonic potentials is discussed in Sect.El Even though no 
(rapid) oscillation is present in the initial data, the solution is rapidly oscillating 
(at frequency 1/e) for any positive time, and focuses at time t = f (Section [21). 
This can be seen by a stationary phase argument applied to the Mehler's formula 
(see [H]), 

(1.3) ^'(^^^) = ^^--^l^^-^ ' )mdy=:U^it)m. 

Perturbations of the harmonic potential by other potentials (sub-quadratic per- 
turbation, see [3D], [13], [H], or super-quadratic perturbation, see [15]) have been 
studied, and in particular the role of these perturbations on the singularities of the 
fundamental solution of the Schrodinger equation. 
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In physics, nonlinear perturbations are considered, for Bosc-Einstcin condensa- 
tion (see |10]). where the harmonic potential is used for its confining properties, 

where N stands for the number of particles and 5 is a coupling constant (in h"^). 
We study precisely the perturbation of (jl.ip with a nonlinear term, 

with (T > 1/n if n = 1, 2, and —r^ < a < if n > 3. We assume that the 

' ' ' n+i 71—2 — 

perturbation of the initial data is small in 

(1.5) E :=iJi(R")nJ^(i/^(M")), 

where the Fourier transform is defined by 

and that / G S. The space S is equipped with the norm 

ll/llE = ||/||L2 + ||V,/|U2 + ||a;/|U2, 
and we assume 11 lis — >0. 

Remark. Initial data with plane oscillations. Let ^0 G M", and introduce 

u%t, x) = u'{t, x-Co sint)e<^-^ ™^*/^ 
Then solves the Schrodinger equation p.4p . with initial data 

u'lt^o^{.f(x)+r%x))e^^. 

Therefore, describing the solution of p.4|) is enough to describe the solution when 
the initial data have plane oscillations. 

We can also prove some results with a focusing critical nonlinearity (2a = 4/n), 

^^g^ iiedtu' + ^e^Au'':^^-^u'~e^\u''\'^/"u', (t, x) e K+ x R", 

We will consider the focusing case only in the one-dimensional situation, and state 
the corresponding results at the end of this introduction. Similar results for the 
multi-dimensional case would be easy to prove. 

The idea of this paper is the following. Initially, the nonlinear term is negligi- 
ble, essentially because the term jw^p'^ is uniformly bounded in suitable Lebesgue 
spaces, therefore it vanishes in the limit e ^ because of the factor e"'^. Mean- 
while, the harmonic potential makes the solution focus near the origin at time 
t = 7r/2, as in the linear case p.ip . When the focusing effects become relevant, 
that is when becomes of order e~"/^, the nonlinear term is no longer negligible. 
On the other hand, if u'^ is localized near x = 0, the term x'^u'^ becomes negligible; 
only the nonlinear term is relevant near the focus. When the nonlinearity is de- 
focusing (Eq. (II. 4p ). the solution u"^ passes through the focus, and the crossing is 
given by the (nonlinear) scattering operator associated to the unsealed Schrodinger 
equation, 

(1.7) + ^AV' = |?M'"V- 
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Since the nonlinearity is defocusing, dispersive effects in (II. 7p are the same as for 
the free equation. Therefore, the solution u'^ leaves the focus along dispersive rays. 
When rays are dispersed, the energy is no longer localized, the nonlinear term 
becomes negligible again and the harmonic potential makes the rule, as before the 
focus (Th. II. 2p . This process can be iterated indefinitely, and each focus crossing 
is described by the scattering operator (Cor. I1.4| l. 

When the nonlinearity is focusing (Eq. ()1.6|) ). and when the mass of / is critical 
(see [23]), the solution blows up near t = 7t/2 (before or after, see Prop. fTTS]) . The 
focusing effects of the harmonic potential first, then of the nonlinear term, cumulate 
and ruin the existence of the solution (Prop. [T3|) . 

In both situations (focusing or defocusing nonlinearity), two distinct regimes 
occur. First, the harmonic potential leads the evolution of the solution, next the 
nonlinear term does so. The two dynamics superpose: they balance each other 
in the case of a defocusing nonlinearity, and cumulate in the case of a focusing 
nonlinearity. The matching of these two regimes occurs in a boundary layer of size 
e around the focus, as in [5], [T5] . 

Formal WKB expansions suggest that with our choice na > 1, the nonlinear 
term is negligible so long as no focusing occurs. We prove that this holds true. 
It would not be so with the choice na — 1; the nonlinear term would be nowhere 
negligible, and we leave out this case. 

On the other hand, we show that the nonlinear term alters the asymptotics of 
the exact solution near and past the (first) focus. More precisely, we prove that the 
caustic crossing is measured by the scattering operator associated to (II. 7p . This 
phenomenon is to be compared with the results of [S], where focusing is caused 
by initial oscillations, and with the results of [3] (see also [T], [2]), where such a 
behavior was first noticed, for the wave equation. In the present case, focusing is 
caused by the oscillations created by the harmonic potential, but the description of 
the phenomena near the focal point is similar. 

The asymptotic state for p.7p we will consider is defined by 

a-8) ^_(x):=^7(.). 

We assume that / G S and that the scattering operator S acts on , with ?/)+ = 
Sip- G S (see Proposition 13. lOp . which is verified in either of the following cases, 

• a > 2-"+V«-+i2«+4 

4n ' 

• II /lis is sufficiently small. 

Assumption 1.1. Our hypotheses are the following: 

• 1 < 71 < 5 and a > 1/2, so that the nonlinearity \z\'^"'z is twice differen- 
tiable. 

• If n — I, we assume moreover that a > 1. 

• If 3 < n < 5, we take a < — ^ . 

• If n < 2, we assume 

2-n + ^n^ + 12?i + 4 

— either a > , 

4ri 

~ or ll/lls < <5 sufficiently small. 

Remark. We could treat the case n > 6 if we replaced the nonlinear term e"'^ \u^ \'^"' 
by F(e"|u^p)u^, with F smooth and 

F{\z\^)<l + \z\^'^. 

Theorem 1.2. Let 2 < r < -^^2- V it- = 1; t<^ke r — oo. Then under Assump- 
tions mi the following asymptotics holds in n L"^ , 
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• IfO<t< n/2, then 



^0 (27r| cost|)"/2 \cost 
• //7r/2 < t < 37r/2, i/ie^ 

g ^ i . . 



^0 (27r|cost|)"/2^+ \cost 
where ip- is defined by m.8\) and -0+ = Sip- . 

Remark. We will prove actually that these asymptotics hold in a stronger sense 
(see Corollary 12.51 Propositions 13.91 and 13. 22]) . 

We can restate this result when time t = 7r/2 is considered as the initial time, in 
place of i = 0. 



Corollary 1.3. Let ip £ J^. Assume that u'^ solves 
(1.9) 



iedtu' + ^e'Au" = ^-u' + {t, x) G 



■"|t=o 



with \\r'^\\s — and ip satisfies the same assumptions as f. Denote ^± = W+^ip, 
£— >o 

where W± are the wave operators (see Provosition \3.1(]\) . Then with r as in Th. \1.2\ 
and under Assumvtions [l~l[ the following asymptotics holds in L'^ (1 , 

• IfO<t<n, then 

^ e^o\ 2TTSmt ^ Vsinty 



• If -TT < t < 0, then 



^ e^o \27rsini / Vsint 



Remark. In 25j, the author considers equations which can be compared to (|1.1|) . 
that is 

' iedtv' + -e^Av' = V{x)v' + U (-) v", 

(1.10) < 1 

I ^ 1 fx\ 

where J7 is a short range potential. The potential V in that case cannot be the 
harmonic potential, for it has to be bounded as well as all its derivatives. In that 
paper, the author proved that under suitable assumptions, the influence of U occurs 
near t = and is localized near the origin, while only the value V{0) of V at the 
origin is relevant in this regime. For times e <C |t| < T*, the situation is different: 
the potential U becomes negligible, while V dictates the propagation. As in our 
paper, the transition between these two regimes is measured by the scattering 
operator associated to U. 

Our assumption na > 1 makes the nonlinear term short range. With our scaling 
for the nonlinearity, this perturbation is relevant only near the focus, where the 
harmonic potential is negligible, while the opposite occurs for e <^ \t\ < tt. In 
this perspective, a new point in our paper (besides the fact that the problem is 
nonlinear) is that we can tell what happens for any time, as stated in the following 
corollary. 
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Corollary 1.4. Suppose the Assumvtion l 1 . 1\ are satisfied. Let fc G N*. Then, with 
r as in Th. the asymptotics of for 7r/2 + (fc — l)7r < t < tt/2 + kir is given, 
in I? n U , by 

^in^—ink^ ____ __. 

uHt,x) - - — ■ -—f^S^i;^ , 

where S'^ denotes the k-th iterate of S (which is well defined under our assumptions 
on f). 

Remark. The phase shift g™^"*"'^?- is present in the hnear case, for Eq. p.ip . 
and is explained in [TT] . On the contrary, the presence of the scattering operator 
S is typically a nonlinear phenomenon, as in [S]. The new point here is that this 
operator is iterated, at each focus crossing. 

Remark. If the nonlinear perturbation was of the form e"*^! jw^p'^^ye^ with cti > 
(72 > (no additional assumption on 172) and nai > 1, the nonlinear term would be 
everywhere negligible, that is, S'^ should be replaced by the identity in Corollarv ll.4l 
This can be seen by an easy adaptation of the proof of Theorem 11.21 This shows 
that the scaling ()1.4p is critical for the nonlinearity to have a leading order influence 
near the singularities {t = Tr/2 + kir). 

We conclude this introduction by stating our result when the nonlinearity is 
focusing (Eq. (|1.6p ). 

Proposition 1.5. Let n = 1 and let R be the unique solution (up to translation 
and sign change) of ~\R" + R = R^ , given by, 

3I/4 

(1.11) Rix) 



cosh {2xV2) 

For t^ G M, define f by f{x) — i?(a;)e*"^^ , and the approximate solution by 

^/ ^ 1 f X \ ,- e ,- 1^1'' 

(1.12) v^{t,x) = =R e -/2+^'*~*e z^-:'*) , 



+et,-t \§+et,~t^ 
Let be the solution of \1.6\) . Then for any A > 0, 

limsup sup \\B'[t){u^{t)-i''{t))\\^2 — * 0, 

£^0 :|-Ae<t<f +et.-Ae A^+oo 

where B^{t) is either of the operators Id, edx or xj e + i(t — n /2)dx. In particular, 
liminf sup \\edxu'^ {t)\\L2 — > +00, 

0<t<7r/2+t.e-Ae A^0+ 

liminf sup ||\/iu'^(Olli°° — +00. 

0<t<7r/2+t.e-Ae A— 0+ 

Remark. The blow up occurs at t = ^ + et^, no matter the sign of t*. This means 
that u'^ can blow up before or after the focus. 

This paper is organized as follows. In Sect. ^ we study the linear equation (|l.ip 
using WKB methods, and introduce two operators ( and H'^) whose role is crucial 
in the nonlinear setting. In Sect. [31 we analyze the nonlinear equation (jl.4p . and 
we prove Th. 11.21 In Sect. HI we prove Prop. fT31 Finally Sect. E] addresses the case 
of anisotropic harmonic potentials. 

Some of the results presented in this paper were announced in [7] . 

Acknowledgment. The idea of studying such equations as p.4p arose from discus- 
sions with L. Miller, after an invitation by P. Gerard and N. Mauser at the Erwin 
Schrodinger Institut in Vienna. This work was partially supported by the ACI grant 
"Equations des ondes : oscillations, dispersion et controle" and the START project 
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"Nonlinear Schrodinger equations" of N. Mauser. The preliminary version of the 
manuscript was improved thanks to remarks made by F. Castella and B. Bidegaray. 

2. WKB EXPANSION FOR THE LINEAR EQUATION 

We seek an approximate solution of the linear equation of the form, 

(2.1) vl^^{t,x)^voit,x)e^^^''^^/'. 

To cancel the term when plugging this approximate solution in (II. ip . the phase 
If must satisfy the eikonal equation, 

(2.2) a,^+i|v.^|2 + M!^0. 

To cancel the term e^, the amplitude vq must satisfy the transport equation, 

(2.3) dtVQ+Wxip-'VxVo + ^VQAip^^O. 

To solve the eikonal equation, one computes the bicharacteristic curves associated 
to the classical Hamiltonian 

p(i,a;,r,0=r + ^ + ^, 

given by 

= 

= 0, 
— 

Therefore, 

x{t) = xq cos t + ^0 sin t, ^{t) — cos t — xq sin t. 
Since no oscillation is present in the initial data, = 0, and the rays of geometric 
optics are sinusoids, 

(2.4) x{t) — Xq cost. 

They all meet at the origin at time t = tt/2, and periodically at time t — tt/2 + kir 
for any keW (Fig. [J). 

Given ^(t) = S/x<f{t), one can solve (|2.2p for < t < 7r/2, by 

Lp{t, x) = — tani, 



and (|2.3p is solved by 

(^^iW^(^)' 

therefore 

I 1^ 

Recall that V{x) — The approximate solution solves 

i <pp|t=o = fix)- 

Denote the remainder := — It solves, 

app ^ 

^\t=o ^ 0- 
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1 1 1 1 1 1 - 

X 

Figure 1. Rays of geometric optics. 

Recall that 

Avo{t,x) = (^) . 

cos-^ t (cost)"/^ \cost/ 

Recall the classical result, 

Lemma 2.1. Assume a function satisfies 

(2.8) iedtw'' + ^e^Aw^ = U{t,x)w'' + S^{t,x), {t,x)eIxR", 

where U is a real-valued potential, I is an interval, and £ Ct{I,L'^). Then the 
following estimate holds for t €z I , 

edt\\w'{t)\\L2 <2\\S'{t)\\L2. 

Applying this lemma, it follows, 

(2.9) edt\\w^t)\\L.<s'^\\Af\\L.. 

COS"^ t 

With the idea of a nonlinear perturbation in mind, it is natural to seek estimates 
in other spaces than L^, in particular Sobolev like spaces. In geometrical optics, it 
is classical to assess e-derivatives to get nonlinear estimates (see for instance |26j). 
This is because e-oscillating solutions are studied. This approach is sharp for multi- 
phase problems, but it contains no geometric information (given by the phase(s)). 
In our case, only one phase is present, and in the nonlinear setting (|1.4p . it remains 
so. In the linear case, this means that controlling in Lebesgue's spaces is 
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equivalent to controlling vq in L^. With Gagliardo-Nirenberg inequalities in mind, 
it is therefore natural to introduce the operator 



(2.10) 



J^it) = -i(cosi)e'"^/^V:r(e""^/^-) = - sint - i costV^. 

s 



Given the dynamics of the harmonic potential, it is also natural to introduce the 
"orthogonal" operator. 



(2.11) 



(t) = X cos t + ie sin tV^ 



When n > 2, we write, for I < j < n, 

(t) = Xj cos t + ie sin tdx^ , 



(2.12) 



Jj (t) — — sin t — i cos tdx^ ■ 



We now state all the properties we will need, including the action on nonlinear 
terms. 

Lemma 2.2. The operators and satisfy the following properties. 
• The commutation relation. 



(2.13) 



H^^(t),^s^, + \s^^^^^ 



= 0. 



• Denote M^{t) = g-^T^t^"*, and Q^it) = e'^i^, then J^{t) and H^t) 
read, for t ^ ^Z, 

(2.14) ,r{t) = ~icoatM^t)\/xM^-t), W{t) = ie saitQ%t)V ^Q^ {-t) 

• The modified Sobolev inequalities. For n — I and t ^ -^Z, 



C 



\\w{t)\\L^ < f \\wmii'\\H^t)wm]!^ 

C olli ' 

For n>2, and 2 <r < define S{r) by 



(2.15) 



S{r) = n ( i - i 



Then for any 2 < r < there exists Cr such that, for t ^ ^Z, 



(2.16) 



\\w{t)\\Lr < 



Cr 



cosi|'5('' 



\wmi-'^^^\\H%t)wm'p. 



|esini|'5('') ' 

For any function F G C^{C,C) satisfying the gauge invariance condition 

3G e C(M+,M), F{z) = zG{\z\^), 
one has, for t ^ ^Z, 



(2.17) 



H^{t)F{w) = d^F{w)H^{t)w ~ dgF{w)H^{t)w, 



J''{t)F{w) = d,F{w)J''{t)w - diF{w)J^{t)w. 
Remarks. 

• Estimates (|2.15p are easy consequences of the conjugation properties (|2.14p . 
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• With the WKB approximation (|2.5p in mind, the | cost|^^/^ term in (|2.15p 
gives optimal time dependence of the L"!^ estimates of the solution of (jl.ip 
away from the focus. This is the main advantage of this operator over all 
the others one could think of (such as Vx in particular). 

• The |esini|~^/^ term in ()2.15p gives optimal L'^ estimates of the solution 
of (jl.ip near the focus (where | smt\ ^ 1). 

• The operator can be considered as the modification of the Galilean 
operator x + itVx, which is very useful in scattering theory (see [H], [TB], 
[T7]). For semi-classical problems where focusing at the origin occurs, it 
was used in 'S] and [B], with the rescaling | + i(i — t*)Va;, where is the 
focusing time. The operator is that operator, transported to the case of 
a harmonic potential. 

• Property (|2.17p states that H'^ and act on nonlinearities satisfying the 
gauge invariance condition like derivatives (Eq. (j2.17p holds for the operator 

• The fact that all these identities, except the first one, hold only for almost 
all i G M is not a problem, since in any case integrations with respect to 
time will be performed. 

• The operators and H"^ are known in quantum mechanics, as Heisenberg 
observables (metaplectic transforms, see [1^], [13]), 



(2.18) r{t) = U^t){~iS/x)U%-t) 

(2.20) H^it) = W[t)xU^{-t) 

(2.21) = U^(t--^){^eVxW[^~t), 



where U'^{t) is the propagator defined by Mehler's formula (|1.3I1 . that is 

The commutation properties (|2.13p are straightforward consequences of the 
conjugation relations (PT^ and ^T^. Identities between ((TTg)) and (P7Tg|) 
on the one hand, (|2.20p and (|2.2ip on the other hand, are due to the geo- 
metric properties of the harmonic oscillator, that rotates the phase space. 
It is easy to check that Sobolev inequalities ()2.15p follow from p.l9p . (I2.20p 
and the estimate 

\\U'it)f\\L^< , I .^m/2 ll/IU" 
(e|smt|)^/^ 

The most remarkable fact is certainly that in the case of the harmonic 
potential, one can estimate the action of these observables of Heisenberg 
on a large class of nonlinearities, through (|2.17p . 

Lemma 12.21 makes it possible to get more precise estimates of the approximation 
given by (the first term of) WKB methods. Denote 

(2.22) ?{:={/€ ij3(R"), such that xf e H'^{W)}. 

Proposition 2.3. Assume f ^ Ji. Then there exists C = C(||/||/f3, ||a;/||//2) such 
that the remainder ~ u^pp satisfies, for < i < 7r/2, 

W - <pp)WIU^+ll^^'(«^ - <pp) WIU^ + W{v' - <pp)(i)IU^ < 

* ds £^ \ 

cos^ s cos^ t I 
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Remark. As mentioned in the introduction (Eq. (|1.3p , the expression of v"^ is given 
exphcitly by an osciUatory integral, and the above result could be proved by sta- 
tionary phase methods. Nevertheless, we do not use this approach, and rather 
present the approach whose spirit is the same as in the nonlinear setting. 

Proof. The first estimate is given by (|2.9p . with C — ||A/||^2. For the second 
estimate, apply J'^{t) to (|2.7p . The commutation property (|2.f 3p yields. 



(2.23) 

One has explicitly, 

~ ^ (c0St)"/2+2^"^^ Vcosj' 

and the same estimate as for the L? case follows, with C = WJWh-^. For the last 
estimate of the proposition, apply H^{t) to (|2.7p . Because of the commutation 
property (|2.13p . the remainder H^{t)w'^ is estimated by the norm of 



H' 



(t) (e-(^-)/^A.o) =(xcos. + ..siniV.) (e^^*-* (^^^^^^A/ (^)) 



cost (cOSt)^^/2+^ \GOSt 

+ ^^tant^— ^i^V.A/(^). 

The norm of the first term is ^^^2 ^ H^A/||^2, and the norm of the second 
term is s ||/|| ^3 . This completes the proof of Proposition l2.3l □ 



From Prop. 12731 WKB methods provide a good approximation of the exact solu- 
tion before focusing. More precisely, the remainder will be small up to a boundary 
layer of size e around t = Tr/2. 

The assumption / S !K is necessary to estimate precisely the validity of WKB 
approximation, but is not really essential. Since the set of such / is dense in S, 
the following lemma shows that this extra regularity can be introduced without 
modifying the asymptotics. 

Lemma 2.4. Assume / e E, and let v"^ be the solution of Then for any 

t > 0, 

Wv^miL^ - II/IIl^; \\J%tXh-2 - IIV/IU.; \\H%t)v-\\L2 = \\xf\\L2. 

Proof. This lemma is a straightforward consequence of Lemma l2.1l and of the com- 
mutation property (|2.13p . □ 

Notice that the L^-norm of v1pp{t) does not depend on time, nor that of J^(i)fapp 
or H^{t)v1pp. We can therefore remove the smoothness assumption of Prop. 

Corollary 2.5. Assume / G S. Then, 

limsup sup WA^it) {v" -vl ){t))\\ — > 0, 

£^0 0<t<f-A£ ^ 

where A"^ (t) is either of the operators Id, J'^ (t) or (t) . 
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3. The nonlinear case 

The proof for asymptotics in the nonUnear setting rehes on Strichartz estimates 
(even though we could do without when n = I). We first recaU how we get them 
in the present case, then prove a general estimate. Then the proof of Th. 11.21 is 
essentially split into three parts: the asymptotics before the focus (0 ^ tt/2 — t), 
the matching between the two regimes (linear and nonlinear) , and the asymptotics 
around the focus {\t — tt/2\ < e). 

3.1. Strichartz inequalities. First, recall the classical definition (see e.g. [5]), 

Definition 3.1. A pair {q,r) is admissible if 2 < r < (resp. 2 < r < oo if 
n — I, 2 < r < oo if n — 2) and 

- = S{r) = n - - - 
q \2 r 

Strichartz estimates provide mixed type estimates (that is, in spaces of the form 
Lj{L]^) with (g,r) admissible) of quantities involving the unitary group e*2^ (see 
EI|,[S|,[22,[lH],[H],[IS],[n])- With the scaling of Eq. ((Lll), the natural unitary 
group to consider is 

(3.1) U^{t) -.^e''^"^. 

Now we can state the Strichartz estimates obtained by a scaling argument from the 
usual ones (with e = 1). The notation r' stands for the Holder conjugate exponent 
of r. 

Proposition 3.2. (Scaled Strichartz inequalities) 

(1) For any admissible pair {q,r), there exists Cr such that 

(3.2) llt^o W«IIl,(R;L-) < Cr\\u\\L^. 

(2) For any admissible pairs (gi,ri) and ((72, 7*2); o,nd any interval I, there 



exists Cri,r2 such that 



1^1 



(3.3) e~^^ 



U^{t- s)F{s)ds 

in{s<t} 



L'!l(/;L'-l) 

The above constants are independent of e and I . 

The proof of this result relies on two properties (see [Sj , or [53] for a more general 
statement): 

• The group C/q is unitary on i^, \\U^{t)\\L2^i^2 = 1. 

• For t 7^ 0, it maps Li(R") into i°°(M"), 

1 



< 



- (e|i|)"/2- 

As a matter of fact, these two estimates also hold for the propagator associated to 
the Schrodinger equation with a harmonic potential (|l.ip . Therefore we can obtain 
similar Strichartz estimates (see [5]). 

If solves (jl.ip . then Mehler's formula yields, for t ^ ttZ (see }12)). 

Therefore: 



The group is unitary on \\U'{t)f\\L2 - 
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For t e] - 7r,0[U]0,7r[, it maps Li(M") into L 

1 



oo /Tipn \ 



(e|sint|)"/2- 

Since for |t| < tt/2, \ smt\ > ^\t\, the proof of Prop. still works when Uq is 
replaced by J7^, provided that only finite time intervals are considered. 

Proposition 3.3. (1) For any admissible pair (g, r), for any finite interval I , 
there exists C'r{I) such that 

(3.4) < Crmu\\L2. 

(2) For any admissible pairs (gi, ri) anrf ((72, ^'2)7 a'^'^ a'^Z/ finite interval I , there 
exists Cri.r2{I) such that 



(3.5) £"1 "2 



U^{t- s)F{s)ds 

lr\{s<t} 



< C',.i,r2(^) II^IIl'2(/;L'-2) 



1 


1 

= - H 


2cr 

1 


i! 


r 


s 


1 


1 


2cr 


9^ 


= - H 
9 





L9i(7;L'-i) 

T/ie above constants are independent of e. 
3.2. A general estimate. We start with an algebraic lemma. 
Lemma 3.4. Let n > 2, and assume < a < There exists q, r, s and k 

— ' 71+2 n—2 ; _? _ — 

satisfying 
(3.6) 



and the additional conditions: 

• The pair {q,r) is admissible, 

• 0< l<S{s) <1. 

If n — I, we take {q,r) ~ (00, 2) and {k,s) = (2cr, cx)). 
Proof. With 6{s) — 1, the first part of p.6p becomes 

and this expression is less than 1 for a < -f^- Still with 5{s) = 1, the second part 
of dSSl) yields 

2 _ n 1 

which lies in ]0, 2[ for < cr < By continuity, these conditions are still 

satisfied for 5{s) close to 1 and (5(s) < 1. □ 

From now on, we assume n > 2 and < a < -^^^2- We state a general estimate 
that can be applied to nonlinear Schrodinger equations with or without harmonic 
potential. Let ll^(t) be a group for which Prop. [331 holds (typically, Uq or in 
our situation). We seek a general estimate for the integral equation, 

u^it) = W{t - ta)ul - ie'""-^ f U'{t ~ s)F''{u'){s)ds 

(3.7) 

-ie-^ / W{t- s)h%s)ds. 

J to 

This equation generalizes the Duhamel formula for Eq. ()1.4p . 

• to the case of the same equation without potential (take Uq in place of C/^), 

• to the case of any initial time and any initial data (ug and to are general). 
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• to the possibility of having a nonhnear term which is not a power, 

• to the possibiUty of having a source term, h'^ . 

Proposition 3.5. Let ti > to, with \ti — to I 5: i"- Assume that there exists a 
constant C independent oft and e such that for to <t<ti, 

(3.8) \\F%u^)m^.' < - — ^-^Ik^WL- , 

(|cost|+e) 

and define 

( t \ 
U' (|co,,t)M.. J ■ 

Then there exist C* independent of e, to and ti such that for any admissible pair 

+ C*e'<'^^'>-i)A%to,t,)\\u'\\mto,t,;L.y 
We will rather use the following corollary, 

Corollary 3.6. Suppose the assumptions of Prop. \3.5\ are satisfied. Assume more- 
over that C*e^'^^^^-^^'^^A'^{tQ,ti) < 1/2, which holds in either of the two cases, 

• < to < ti < ^ — Ae, with A > Aq sufficiently large, 

• to,ti E — Ae, ^ + Ae], with < rj sufficiently small. 

Then 

(3.10) ||^^1U~(to,ti;L^) <C||<|U.+C,,,e-i-^||/^^|L,-(,„,,^^i.-). 

Proof of Proposition 1 3. 5[ Apply Strichartz inequalities p.4p and (|3.5p to p.7p with 
qi = q, ri ~ r, and q2 ~ q, r2 — r for the term with F^(tt^), q2 = q, r2 — r for the 
term with h^ , it yields 

\\u'\\mtoM;L^) <Ce-i/2|l^^|U.+C,,,e-'-i-'|l/i^|L,.(,„,,^;i..) 

Then estimate the space norm of the last term by (13.81) and apply Holder inequality 
in time, thanks to (|3.6p . it yields p.9p . □ 
Proof of Corollary \3.6l The additional assumption implies that the last term in 
p.9p can be "absorbed" by the left-hand side, up to doubling the constants, 

(3.11) <Ce-i/i||ug|U.+Ce-'-i-'||/i^|L,.(,„,,^^^,.). 

Now apply Strichartz inequalities p.4p and (|3.5p to p.7p again, but with gi = oo, 
ri = 2, and 172 — q, ?'2 — L for the term with _F^(u^), q2 = q, r2 ~ r for the term 
with ft,^. It yields 



Like before, 

no"— 1 



ll^'("')llLl'(to,t.;L.') <^^^^'"^''-^"*^^'(^0,tl)h1U£(to,*,;L^) 
<Ceih^|U^(i„,,,,i.), 

and the corollary follows from (|3.1ip . □ 
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3.3. Existence results. Local existence in S steins from the well-known case of 
the nonlinear Schrodinger equation (|1.7p . once we noticed that the operators 
and J"^ are the exact substitutes for the usual operators eV and | + z(t — ■|)V, by 
Lemma [221 Duhamel's formula for (11.41) writes 



(3.12) u'{t) = U'{t){f + r')-ie'"'-^f {t ~ s){\u'f''u'){s)ds. 

Ja 

Replacing If^ with Uq would yield the Duhamel's formula for the same equation 
with no harmonic potential. From the above remark (the essential point is that H"^ 
and commute with C/^) and the fact that the same Strichartz inequalities hold 
for and Uq when time is bounded, local existence is actually a byproduct of the 
existence theory for (|1.4p (which relies essentially on the results of Sect. 13.21 see 
[5^ . [S], [IS], [Uj). For (go, To) admissible, introduce the spaces 

= K e C(/,S),V(q,r) admissible, 7^^^t^ J'^m'^ e Ll^{I,Ll)}. 

Proposition 3.7. Fix e e]0,l], and let f,r^ e S. There exists t^ > such that 
([j.^[ ) /las a unique solution £ ^2cr+2(0j ^^)- Moreover, this solution belongs to 
y^(0, f^). The same result holds for Eq. il.b]) and for any initial time. 

We can take t^ = +00 when the nonlinearity is defocusing (Eq. (II. 4^ ). thanks to 
the conservations of mass and energy, 

(3.13) \\u^t)\\L2 = \\u^0)\\L2=O{l), 

E%t):^-\\eV,u%t)\\l2+ / V{x)\u%t,x)\'dx+ —-\\u%t)\\f+^2 

(3.14) ^ jm" ct + i 

^E'{0) = 0(1). 

The conservation of energy provides an a priori estimate for H^u'^ and J^u^ thanks 
to the identity, 

(3.15) Vt,a;, \H' {t)u^ {t,x)\^ + {t)u'' {t,x)\'^ = \x\'^\u''{t,x)\^ + \eV ^u^{t,x)\'^ . 

Proposition 3.8. Fix e g]0, 1] and let f^r^ G S. Then |i.^[ ) /las a unique solution 
y,e g y£(K) t/iere exists C such that for any < > and any e g]0, 1], 

(3.16) \\u'{t)\\L^ + \\£^xu'{t)\\Li + \\xu'{t,x)\\Li < C. 

3.4. Propagation before the focus. Before the focus, we take as an approximate 
solution the solution of the linear problem, that is, defined by p.ip . 

Notice that from Prop. 12.31 we know the asymptotic behavior of before the 
focus. We prove that in the very same region, is a good approximation of the 
nonlinear problem. 

Proposition 3.9. Assume f,r'^ e S. Then 

limsup sup \\A''[t)[u''[t) - v^{t))\\^2 — > 0, 

£^0 0<t<f-Ae A^+00 

where A'^{t) is either of the operators Id, J'^{t) or H^{t). 
Proof. Define the remainder w"^ = — v'^ . It solves 

iedtw' + ^e^Aw' = Vix)w' + e"''\u'\^''u' , 
From Duhamel's principle, this writes, 

(3.17) w'{t) = U'{ty - le'""^ f U'it - s) (s)ds. 

-'0 



SEMI-CLASSICAL NLS WITH HARMONIC POTENTIAL 



15 



Since v'^ solves the linear equation (jl.ip . so does J^{t)v^ , and 

\\v%m^.^\\fh., iiJ^(iKiu. = iiv/iu.. 

From Sobolev inequality (|2.16p . 



C 



l-'5(l)||Y7f||'5(£) 



Therefore there exists Co such that 

(3.18) ||T;"(t)|L, < , ^° , , . 
^ ^ " ^ - |cost|*(^) 

From Sobolev inequality, for e sufficiently small, ||u'^(0)||^3^ < Cq. From Prop, 
for fixed e > 0, m*^ e C(M, E), and the same obviously holds for w^. Therefore, there 
exists > such that 

(3.19) < , ^° , , , 
^ ^ II ^ ^iii- - Icosil-'t^) 

for any t £ [0, i^]. So long as p.l9p holds, we have 

2Co 



\u%t)\\^,< 



urn' 



and we can apply Prop.EH Indeed, take It"^ = ^7^ h'' = e""" {u" v"" and F^(w^) = 
From Holder inequality and the above estimate, 

(|cosi|) 

Assume ([XTg]) holds for < i < T. If < t < T < f - Ae, then e < cos t, and the 
above estimate shows that F'^ satisfies assumption p.Sp . 

From Cor. 13.61 if A is sufficiently large, then forO<t<r<^ — Ae, and for 
any {q, r) admissible, 

\\w^\\L^io.T;L^) < Cllr^lk^ + Ce-'^-i-^ II I"1"^^1L.'(o.T;L^') ■ 
Taking {q,r) = {q,r) yields, from Holder inequality, 

l|l'^'l'"^iLl'(0,T;L^') ^ h"\\lHO,T-M)\\^"\\LHO.W- 

The first term of the right-hand side is estimated through p.lSp and (|3.19p . The 
last term is estimated the same way, for p.lSp still holds when replacing s with r. 
Therefore, 

l|KP'^-1L.'(0,T;L.')<- ^ 



and 



/ \ 'na—i 

(3.20) <C|lr^|U.+c(^^^j 

Now apply the operator to (|3.17p . Since and t/^ commute, it yields, 
r{t)w' = U'it)J''{Oy - ie"^-i / U'{t - s)J'{s) {lu'l^^u") {s)ds. 
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Because acts on this nonlinear like a derivative, we have an equation which 
is very similar to (|3.17p . with w"^ replaced by J^w^ and replaced by — iVr^. 
Therefore the same computation as above yields 

(3.21) II JW|Uoo(o,T;L^) < C||Vr^||i2 + C 



2 ^ 



Combining ((X^ and yields, along with ((TTC)) . 



yte[Q,T], \\w^t)\\^^< 



c 



^ nZTTITIT II' 11^^ 



Therefore, choosing e sufficiently small and A sufficiently large, we deduce that we 
can take T = f - Ae. This yields Prop. for = Id and The case A" = 
is now straightforward. □ 

3.5. Matching linear and nonlinear regimes. When time approaches 7r/2, the 
nonlinear term cannot be neglected. On the other hand, since the solution tends 
to concentrate at the origin, the potential becomes negligible. It is then natural to 
seek an approximate solution that solves 

iedtv' + ^e^AS^ = £"'"|rp'^r. 

The question that arises naturally is, how can we match and v^7 With the 
results of [F in mind, we can expect that is exactly a concentrating profile, 

(3.22) r(t,a:) = ^^(^- 

The function -0 must be defined to match the solution , or one of its approxi- 
mations v'^ or v^pp, when t = tt/2 — Ae, for A sufficiently large. Notice that this 
problem was already encountered by H. Bahouri and P. Gerard in [3] (see also [T], 
I. Gallagher and P. Gerard [H]). We prove that for A > sufficiently large, the 
propagation for 7r/2 — Ae < i < 7r/2 + Ae is described by w^. 

Write tl = 7r/2 — Ae, and assume from now on that A > 1. For large A, Prop. [23 
and 13.91 imply 

u^itl.) ^ v^itl.) ^ vlppitl.) . j-^f (-^) e-^. 
For Ae close to zero, the following approximation is expected, 

1 f X \ \ / X \ - 

(sin(Ae))"/2'' \sm{Ae)J (Ae)"/2-' \AeJ 

We prove that this approximation is correct in Lemma 13.131 below. From (|3.22p , 
this should also be close to 



^n/2 

Recall the classical result. 

Proposition 3.10. d^U], Theorem 1.1; [5], Theorem 4.2) Assume ip- € S and 

n+2 ^ n-2 



—To < tr < — i-n if n>2, a > 1 ifn=\. Denote 



, , 2 - n + Vn2+T2nT4 

Tn ■ 

If cr > aQ{n) or if Wip-Ws is sufficiently small, then 

• There exists a unique ijj G C(M(,S) solution of |j. 7| ), such that 

lim IIV'- - t/o(-i)V'(^)l|s = 0, whereUait) ^ e'^"^. 
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• There exists a unique "0+ S S such that 

lim U'+ - Uo{-t)^Kt)h = 0. 

t — *+oo 

Recall that the asymptotic state ■(/'_ was defined in introduction by, 

1 



V'- 



r/, 



(2i7r)»/2- 

and the approximate solution (near t = 7r/2) is given by 
We prove, 

Proposition 3.11. Assume f^r'^ € S. Take tp^ defined by il.S]) . Then 



lim sup 



Ae, . 



^(C/o(-A)V.-)(-) 



L2 



A^+oo 



0, 



and the same holds when applying either of the operators eV^ or - — iAsVx to the 
considered functions. 



Proof. From Cor. 12.51 (from which ~ "^app) Prop. 13.91 (from which 



(3.23) lim sup 

£—0 



u%tlx) 



-appy 

1 „ f X \ liL 



g '2etan(Ae) 



L2 



A^+oo 



0, 



(sin(Ae))"/2 \sin(Ae)^ 

and the same result holds when applying either of the operators J^{tl) or H^{tl). 
Notice that applying J'^{t%) or H'^{t%) is not so different from applying eVa, or 
I — iAeWx, for when Ae goes to zero, 

J%tl)^--iAeVx , H%tl)r^^sVx■ 
e 

Recall that = 7r/2 - Ae. 

Lemma 3.12. Let a^(t^, .) G E be a family of functions such that there exists C* 
independent of e g]0, 1] such that, 

(3.24) \\xa%tl,x)\\L2 + \\eVxa%tl,x)\\^, < C,. 

Then for any A > 1 , 



lim sup 



(j^'(t:)- J+iAeV,) a'itl) 



L2 



limsup||(i7''(t^) - ieWx)a''{tl)\\L2 = 0. 

£^0 



In particular, we can take ~ u'^ or — ffpp- 

Remark. Lemma 13.121 has a simple geometric interpretation. Near the focus, rays 
of geometric optics, given by ()2.4p . are straightened (Fig. [5]). Thus in the neighbor- 
hood of t = 7r/2, rays are almost straight lines, that is, the geometry is nearly the 
same as in [5]. In that case, with the natural scaling p.22p . the "good" operators 
are eVx and - + i{t ~ 'k/2)Vx- 



Proof of Lemma[3JE Fix A > 1. 

(^J'itl) - ^ +iAeV,) a'{tl,x) = ((cos(Ae) - 1)^ - z(sin(Ae) - Ae)V,) a'{tl,x) 
Taking the norm yields, 
(j'{tl)--+iAsS/x)a'{tl) <C{Aef -a'{tl,x) ^C(Aef\\Vxa'{tl)\\ 
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Figure 2. Rays of geometric optics are straightened near t — ^. 



The assumption (13.24p (which is a consequence of (|3.16p for , and straightforward 
for Wapp) implies 



L2 



which proves the first part of the lemma. Similarly, 

< C{Ae)+C{Ae f. 

This completes the proof of the lemma. □ 



Now we prove that in p.23p . we can replace sin(Ae) and tan(Ae) with Ae up to 
a small error term. Denote 

Lemma 3.13. Assume / G S. For any A > 1, 

limsup||«pp - w^pp)(i^)||^, lim sup II eV:r«pp -w^pp)(t*) 11^,2 = 



lim sup 

£-►0 



-*AeV,) «pp-wipp)(t^) 



0. 



Proof of Lemma \3.1S[ Write A = Ae. For fixed A, A is a small parameter when e 
goes to zero, and 

^(sinA)"/2^ VsinA/ A"/2-'vA/y 
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Taking the norm yields, 



inA)"/2-^ (sinA) A"/2-^(a) 



(sin A) 



< 



(sin A)"/2 A"/2 
/ 



< 





1 


+ 


A»/2 


+ 


/(^) 


f sinA^ 


I A J 


+ 





sin A 



sin A 
A 



L2 



L2 



L2 



e 

n/2 



2e V ^ 



1 



L2 



l/ll 



sin A 



-/(•) 



2e V tJi^ri A 



1 



L2 



The first term of the right-hand side clearly goes to zero with A. So does the second 
one: if / G C^(M), it is 0{X^), and by density, it is o(l) when A goes to zero for 
any f E L^. Recalling that A = Ae, we have 



A2 



1 



e V tan A 



-A = A 



A 



tan A 



- 1 



Thus, for any fixed A > 1, this term goes to zero when e goes to zero. Therefore, 
from dominated convergence, for any fixed A > 1, 



lim sup 



1 



0. 



L2 



Computations for ||£V^«pp - wlpp)(t^)||^2 and ||(f - iAe\7^) «pp - wlpp)(<^)||^2 
are similar and essentially involve one more derivative or one more momentum. In- 
deed, vlpp and vlpp behave well with respect to the operators e\7x and | — iAeV^:, 
thus we can use the same density argument as above. □ 

The next step to prove Prop. 13.111 consists in comparing v^pp and the rescaled 
free evolution of the asymptotic state ip- . 

Lemma 3.14. Assume / G S. The following limits hold, uniformly with respect to 
eG]0,l], 



lim 

A—y+rx 



-ri/2 



([/o(~A)V'- 



L2 



= lim 

A^+oo 

= lim 

A^+oo 



-app(tD-^(C/0(-A)^-)(-) 

)(^app(^*)-^(t^"(-A)V^-)(; 



= 0. 



Proof of Lemma \3.14\ From the Fourier Inversion Formula, we have 

1 1 



_ ■ Jil_ 



(27r)" (Ae)"/2 



e-'^ / eTKf{y)dy. 
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On the other hand, the expression of the free Schrodinger group Uq imphes, along 
with definition ()1.8p . 



lj([/„(~A)^_)(^) 



27rAe 
1 



ri/2 



1 



(27r)" (Ae)"/2 
Thus the remainder we have to assess writes 
1 1 w 



■y Iy , , s , 



^'^ f{v)dy. 



(27r)" (A£:)»/2 



which is also, 



(( 



1 



1 - e*5A / 



^ f{y)dy, 



(^) 



(Ae)"/2 

The lemma then follows from the strong convergence in L^, e^^^ 
Lemmas [XT^ [XHl and [XTH imply Prop. [TTT] 



5^0 



□ 
□ 



3.6. Description of the solution near the focus. Propositions 13.1(11 and [TTT] 

imply that 



lim sup 



- - Ae, • - ( - - Ae, 



L2 



lim sup 
lim sup 



eV, (u^ (I - Ae, •) - (I - Ae, 



a; 

zAeV 



L2 



A^+oc- 



0. 



This means that taking A large enough, and e small enough, the difference — 
is small at time = 7r/2 — Ae, which is the "initial" time in the boundary layer 
where nonlinear effects take place (and where the potential is negligible). Since the 
role of is negligible, we first assume = Q. 

Proposition 3.15. Assume f £ IK, and that the nonlinearity is , that is, a > 
1/2, which is possible only if n < 5. Then the difference u'^ — v"^ is small around 
the focus. 

limsup sup - — > 0, 

e^O f-Ae<t<f+Ae A^+oo 

where (t) is either of the operators Id, J^{t) or H^{t). 

Remark. The assumption cr > ^ is needed to prove Lemma 13.171 below. It seems 
purely technical, and one expects Lemma 13.171 to hold without this assumption. If 
n = 2, the nonlinearity is automatically thanks to the assumption a > If 
n = 3, then we have to restrict our study to the case ^ < cr < 2. In particular, the 
value CT = 1, which corresponds to a cubic nonlinearity, is accepted. 

Proof. Propositions 13.101 and [XTT] imply that 



lim sup 

£^0 



^ A 



A 



lim sup 
lim sup 



ev(.^(|-Ae,.)-rg-Ae,. 



L2 A^+oo 
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Define the remainder w'^ ~ u"^ — , and keep tlie notation t% = 7r/2 — Ae. From 
Prop. [TTT| 

\imsn^\\B'{tlW{m^, 0, 

where is either of the operators Id, - + — 7r/2)V or eV. From Lemma l3.12l 

this implies 

Yiu,sM\A'{tl)w'(tl)\\^, ^ 0, 

where A'^{t) is either of the operators Id, J'^{t) or H^{t). 
From the conservation of energy (|3.14p . we have 

\\eWu'{t)\\L2<C. 

From the conservation of energy for (|1.7p . 

I (^iiv^(i)iii. + ^wmwT^t^^) = 0, 

we have 

\\sVv%t)\\L2 < C. 

Therefore, since = u"^ — v'^ , 

\\eWw'{t)U2<C. 

From Sobolev inequahty, 

and there exists Ci such that for any t e M, 
(3.25) l\^^t)U.<^. 

This estimate will be useful for \t — 7r/2| < Age, where Aq is given by Cor. 13.61 
For \t — n/2\ > Age, sharper estimates are provided by J^, along with Sobolev 
inequality (|2.16p . 

The first step of the proof consists in showing that the harmonic potential can 
be truncated near the origin without altering the asymptotics. Let x G C^(M") be 
a cut-off function, with 

suppx C B{0, 2), < X < 1 and Vx G B{0, 1), xi^) = 1- 
For R> 0, define 

URit,x) =x[^) u%t,x). 

Lemma 3.16. Assume f £ "K, a > ^ and take i? = e". Then for any < a < 1, 
limsup sup \\A'{t){u''{t) -u]^{t))\\L2 — > 0, 

E->0 |_Ae<t<f+A£ ' A^+oo 

where A'^{t) is either of the operators Id, J'^(t) or H^{t). 
Proof of Lemma \3.1b\ The function ujj satisfies, 

therefore the difference w|j := — u'j^ solves, 

m + 'r'A - M!) . ,...|„.|..„. _ ilv, (£) v.. - A, (£; 

From Lemma |2.H and because the term £"'^|u^p°' can be considered as a real 
potential, 
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which imphes, from (|3.13p and (|3.16p . 

Integrating this inequahty on — Ae, ^ + Ae] gives 

Ikl^WIU^ < \\w%{itI2 - A£)|U. + CA| + CK 



sup 

-A£<t<^+A£ 



Taking R — e" with < a < 1 yields, 

limsup sup < hmsup ||ii;|j.(7r/2 — Ae)|lL2. 



e^O f-Ae<t<f+A£ 



Now since ip- G i^, < a < 1 implies, along with the dominated convergence 
theorem, 



£^0 



-0. 



From Prop, [nun the first part of Lemma r3.16l (with = Id) follows. 



To estimate J'^u'lj, notice that 



r{t)wUt, x)^[l-x (I)) J'{t)u%t, x) + Vx (I) u%t, x), 
and for f - Ae < < < f + Ae, 



cost. 



R 



vxy.^(^,.) 



L2 



R R 



Therefore to prove Lemma 13.161 when A"^ — J"^ , it is enough to prove. 



lim sup sup 



= 0. 



L2 



The function J^{t)u^ satisfies, from the commutation property (j2.13p . 



(3.26) 



ledt + -e^A - ^ ) .r{t)u' = e-r{t) (| 



Notice that from Prop. [XH and p.l6p . Sobolev inequality implies that there exists 
C = C(A) such that for any t e [0, tt/2 + Ae], 



(3.27) 



C 



(|cosi| + ey(^y 



At this stage, C might depend on A (even though we will know it does not, after- 
ward). Therefore, Cor. 13.61 applied to p.26p a finite number of times to cover the 
interval — AqC, ^ + Ae], implies that for any A > Aq, J^{t)u^ is bounded in 
for t E [0,7r/2 + Ae]. Next, commuting the cut-off function x with (|3.26p yields. 



^sd, + ;-e'A-^]{l-x(^))J%t)u^^e-^{l-x{^))r{t) 



\R 



( ^ 



\R 



,£ \2<7 ^,e\ 
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From Cor. [XI] and (P^Tf]) . if we denote ii = 7r/2 - As, = 7r/2 + Ae, we have 



< 



L°°(ti,t5_;L2) 



< 



1-X 



L2 



We can conclude with the foUowing lemma, whose proof is postponed to Sect. 13.71 

Lemma 3.17. Assume f £ 'K and a > 1/2. Let A > 1. There exists C = C(A) 
such that for any t Cz [7r/2 — Ae, 7r/2 + Ae], 



\\eWr{t)i 



leWH^ity 



< C. 



L2 -r ||C V J2 yLJU. 11^2 

This completes the proof of Lemma 13.161 the computations with being sim- 
ilar. □ 

To prove Prop. 13.151 we now have to compare and the truncated exact solution 

Lemma 3.18. Assume f G 'K and take R — . Then for any < a < 1, 
limsup sup \\A^{t){u'^{t) - i'{t))\\L2 — > 0, 

£—0 f-A£<t<f+Ae A^+oo 

where A'^{t) is either of the operators Id, J'^(t) or H^{t). 



Proof of Lemma \3.18l Denote wf 



— u'^. Recall that mJj solves 



iedt + ie^A 



2 



and notice that with our choice for the cut-off function 

^©=^(2]?)^(S)' 



therefore 



1 



n<j I ^ ,£ 1 2a"^ ,e 



2R 



where 



The remainder solves 



1 



v«(x).,(i)M; 



(3.28) 



iedt + ttE^A w% ^Vr{x)u% + e 



Apply Prop. [231 with now U'' = f7|f , F'^ = and h^ = h\ + hl, where 



2R 



and 
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This yields, for 7r/2 - Ae < io < < 7r/2 + Ae, 



Ce 



<Ce -^WR{to)\\ 



Ce 



-i-i 



C 



h -tr 



from Holder inequality. Taking 



C 



2(7 /fc 



<1 — to 



K^flllL9(to,ti;Li^) 



"R\\L3.(to,ti;Ll^] 



2a I k 



< 



we have 



+ ll^^l?,llLa.(fo,ti:Li:) ■ 



Repeating this manipulation a finite number of times covers the whole interval 
t G [7r/2 — Ape, 7r/2 + Aqe]. Doing this, we get a possibly large, but finite, constant, 
which can be seen as the analogue of the exponential term in Gronwall lemma. For 
AqE < |i ~ 7i'/2| < Ae, we use time decay estimates provided by J^; bearing the 
comparison with Gronwall lemma in mind, this means that the operator provides 
some integrability for Age < |i — 7r/2|, which is stated in (|3.8|) . and implies the first 
condition in Cor. 13.61 This integrability is needed to get a bound independent of 
A > Aq. When A'' = , from ((TT^ . 



iex sint, 



and H^u^j^ satisfies, 



R 



It follows that the remainder H^w'j^ satisfies. 



Vr{x)H^u% + ie(sint)VVR,(x)ufj 
|^Vx(|)i/^a..^ + zesin.^Ax(|)v.^ 

2 

Ax(^) H^u" + iesint^VAx(^) ■ 



+ ixex (^J^ (sini)u^ + i£(sint)VVR(x)u|. 
+ VHix)x (I) i/^z^^ + ^^Vnix)isint)Vx (|) 



H^Wu'^ + iesini 



2i?2 

2 
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We can estimate the term in H'^dxu'^ because we can estimate dxH^u"^ fLemma l3.17p 
and the foUowing holds, 

[H^{t)y] = -cost = 0(e) for 7r/2 - Ae < t < 7r/2 + Ae. 

The proof then proceeds as above. □ 

Lemmas 13.161 and 13.181 clearly imply Prop. 13.151 □ 

The assumption / G J{ turns out to be unnecessary. Indeed, we can use a density 
argument for , and approach / g S by functions in J{ up to a small error in the 
norms that are considered in Prop. 13.151 this stems from global well-posedness of 
p.7p (see e.g. [5], [U]). We can mimic the proof of this result for u'^, thanks to 
and H". 

Proposition 3.19. Provosition \3 . 1 5\ still holds if we assume / € S and ^ 0. 
3.7. Proof of Lemma 13.171 We first use the following remark. 
Lemma 3.20. Assume a function u'^{x) satisfies 

IK-e^A + 1x12)^,^11^, + \\eWu'U2 + \\xu'\\l2 < C, 
where C_ does not depend on e. Then 

\\e^Au'\\L2 + \\\x\\'\\L2 <C. 



Now the idea is to differentiate ()1.4p with respect to time. This is classical for 
the case of the nonhnear Schrodinger equation (|1.7p . see e.g. [5], Sect. 5.2. Thanks 
to the above lemma, we can adapt the mentioned results to prove the following 
proposition. 

Proposition 3.21. Assume f e'K. Let A > 1. Then 

e C (o, I + Ae; n T{H^)^ n (o, | + Ae; 

and there exists C = C'(A) independent of e such that 

sup \\edtu^{t)\\L2 + sup \\e^l^u^[t)\\L2 + sup \\\x\^u^{t)\\L2 < C. 

0<t<f+Ae 0<t<f+Ae 0<t<f+A£ 

Idea of the proof. As in , the idea of the proof consists in differentiating the equa- 
tion satisfied by with respect to time, and estimate edtu'^ . Since the harmonic 
potential commutes with the time derivative, one can mimic the proof given in [5], 
Sect. 5.2. When there is no potential, like in (jl.7p . the control of the nonlinear 
term and the time derivative give some control on e^Aw^. In our case, this controls 
e^Au^ — From the above lemma, this means that we can estimate each of 

these two terms. □ 

Notice that the following algebraic identity holds point- wise, for any j, fc, 
(3.29) 

\xkeJ^itX\^ + \s^dkJ^{tX\'' + \xkH]{tX\'' + \edkH]{tX\^ = 

= \xJXkU^t)\^ + \exkd,u^t)\'' + \eS,kU^t)+exjdkU^t)\^ + \e^d%u^t)\\ 

where Sjk stands for the Kronecker symbol. From Prop. 13.211 the right-hand side 
is bounded in L]., uniformly for < t < -I + Ae. This implies the boundedness of 
eVH^{t)u'^ stated in Lemma [3. 17) and even a little more, that is, 

(3.30) Vte 0, l+Ae , \\e\I H^{t)u^\\L2 < C. 

At this stage, we have not assumed that the nonlinearity was twice differentiable. 
On the other hand, we just have 

\\eyj^{t)u'\\L^ < -. 

£ 
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The idea is that for this term, (|3.29p is far from giving a sharp estimate. Indeed, for 
|i — 7r/2| = 0{e), we guess that the main contribution of u'^ hes in |a;| — 0{e) (semi- 
classical Schrodinger equations are morally hyperbolic). This is precisely what we 
have to prove. With the additional remark that near t — tt/2, one can replace 
with eV up to a small error term, this suggests that the leading order term of the 
left-hand side is \£dxkHj{t)u^\'^, and the leading order term of the right-hand side 
is |e^9^^2,^tt^(i)p. Thus there is nothing more to hope from this identity. 

This in fact must not be surprising. The only additional estimates we obtained 
are those stated in Prop. 13.211 

sup \\edtu''{t)\\L2 + sup \\e^Au'{t)\\L2 + sup \\\x\'^u''{t)\\L2 < C. 

0<t<f+Ae 0<t<f+Ae 0<t<f+A£ 

The boundedness of the first two terms means that is e-oscillating, and the 
boundedness of the last term means that the solution remains confined. This is 
due to the fact that we work with an unbounded potential, but not to the fact that 
we consider the harmonic potential in particular. Therefore, there is no precise 
geometric information in this estimate. As a matter of fact, away from the focus, 
this kind of information is given by the operator J^. 

We assume that the nonlinearity is twice differentiable. Recall that for I < j < n, 
Jj{t)u'^ satisfies 

2' 2 



iedt + -e^A 



Differentiating this equation with respect to Xk yields 

(3.31) (^ledt + ie'A - edkJ]it)u' = e^^^+'d^J^it) {\u'\-"'u')+exkJ] it)u' . 

The last term comes from the commutation of the harmonic potential with dx,. ■ 
From (|2.17p . the following point- wise estimate holds, 

(3.32) \edkJ^it) ilu'f'u')] < \u'f''-'\edku'\.\J]{t)u'\ + \u'f''\edkJ]{t)u'\. 

The idea is that the last term is well prepared to apply Gronwall lemma. For the 
first term of the left-hand side, we have to work a little more. Apply Cor. 13.61 to 
(P3T|) . with 

\F'{edkJ]it)u')\ < \u'\^''\edkJ]{tX\, 
\h'\ < e"nw1'""'k9feU^M4{iKI + \exkJ^{t)u'\. 
We already know that F'^ satisfies (|3.8p . Cor. 13.61 vields. 

\\£dkJ^it)u'\\^^ < \\edk{x,f)\\L2+C\\xkJ^{tX\\ 

(3.33) 

+ Cs"--^-i|||.1--sa.uV|(t).i^,,(^^^^^^^,,,. 
For fixed t, Holder inequality yields. 



] ^ ' ML. 

with 



2a - 1 1 2cr 



ai 02 s 

We can take for instance ai = 02 = s. This implies, along with (|2.16p . since 
\\£dkH^{t)u'^\\]^2 is uniformly bounded 
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Now apply Holder inequality in time, with 

1 _ 2(7 1 1 

c/ k CO q 

This yields 

\\\u'\^'^-'edkU-J^{tX\\^,, , < A%to,h)\\J]{tX\\mtoM;L^) 

(3 34) 



where A'^ is defined in Prop. 13.51 We also know that 

\\J]it)u^\mt,,tuL^) < Ce-i/i, 

therefore (|3.33p yields, 
But from Lemma [3^ 

ncr - 1 - I = 2cr (^S{s) - , 

and we find the same quantity as in Prop. 1X51 that is s'^"^^'"-''^^^ A"^ (to,ti). With 
the remarks that 

[edk, Jj{t)] = 5jk sint, 

and S{a2) < 1, we have also, 

(3.35) 2 



Now it is natural to study XkJj{t)u^ . It satisfies, 

= e'^'xkJ^it) {\u'\'"'u') +e^dkJ]{t)i 



The same computation as above, minus the three terms estimate which is not 
needed here, yields 

(3.36) 

+ Ce"" ^A'{to,ti)\\xkJ^it)u'\\L^(^t„^tuL^). 

Summing p.35p and (|3.36p over j and k, Lemma 13.171 follows from the Gronwall 
lemma. □ 

3.8. Past the first focus. After the first focus, we can proceed like before the 
focus, and iterate this process. Notice that if n > 3, then i > (Jo{n), and we always 
have a > ao{n). Next, we can prove the analogous of Prop. [3T9l using Prop. [STTT] 
and Corollarv l2.5l 

Proposition 3.22. The following asymptotics holds for 7r/2 < t < n, 
limsup sup \\A''{t){u'' -vl){t)\\L2 — > 0, 

e^O 7r/2+Ae<i<7r A^+oo 

where A'^{t) is either of the operators Id, J'^(t) or H^{t). 
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Finally, ufpp i approximates vf like in Corollary 12.51 Corollary 12. 5[ Prop. 13.91 
EH and [3121] then imply Theorem O 

When t — IT, the problem is almost the same as at time t = 0. The initial data 
/ is replaced by /i , and 

Therefore, Theorem ll.2l can be iterated, which yields Corollarv ll.4|, because of the 
property, 

ye e M, VV'- e S, 5'(e*V-) = e'^S'(i/'-). 
4. When the nonlinearity is focusing 

In this section, we assume n — 1 for simplicity. The first remark to guess the 
result of Prop. [T751 is that in the proof of Prop. [?75l the sign of the nonlinearity in 
unimportant. One needs local existence results to start the "so long" argument, 
and general estimates on the nonlinear term that do not involve its sign. Therefore 
Prop. [311] still holds when is the solution of (|1.6p . 

Next, assume for a moment that the matching argument can be used as in 
Prop. [5TTTl and that afterward, the harmonic potential can be neglected because of 
concentration. The behavior of should then be the same as the solution of 

Resuming the scaling p.22|) , we have to understand the behavior of the solution of 
the same equation with e = 1. It is well known (see [8]) that for small initial data, 
the solution exists globally. The critical mass is the L^-norm of the ground state 
R defined in Prop. 11.51 Recall what happens in this critical case. 

Theorem 4.1. (]24] . case n = 1) Let e E, with ||</?||l2 = ||i?|j2^2. Let ip be the 
solution of the initial value problem, 

Assume that ip blows up at time t — t^. Then there exist 0, oj, ^q, xi G M such that 
for t < t^, 



The second important remark is that such profiles as in (j4.ip are dispersed when 
t goes to — oo. U uj = 1, xi = £,0 — 0, then 



(4.2) 



Uo{t, - t)ij{t) - -^i? 

V2i7r 



— > 0. 

, t — > — QO 



From the uniqueness in the first part of Prop. 13.101 if i/j solves the critical nonlin- 
ear Schrodinger equation and behaves asymptotically when t —oo like the free 
evolution of 

^- -i=C/o(-t*)i?, 
V2i7r 

then i/j is given by (|4.ip with oj ~ 1 and xi = £,o — 0. Back to the scaling (|3.22p . 
this yields the definitions f{x) = _R(x)e*^l^l (from (jl.Sp and the definition of 
and (fm| . 

Now sketch the proof of Prop. 11.51 As we noticed. Prop. 13.91 describes the 
behavior of up to i = 7r/2 — Ae for large A. What prevents us from mimicking 
the proof of Prop. 13.111 .'' The limit (|3.23p still holds, as well as Lemmas 13.131 and 
13.141 However, one cannot apply Lemma 13.121 so easily to and v^pp for estimate 
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p.l6p is not true when the nonhnearity is focusing. On the other hand, p.l6p is 
true up to time t = 7r/2 — Ae for large A, from Prop. and the algebraic identity 
(|3l31) . Therefore Prop. [SUlstiU holds. 

Finally, one can adapt Prop. lSTSl bv replacing the time interval [tt/2 — Ae, tt/2 + 
Ae] by [tt/2 — Ae, n/2 + t^e — Ae] , for any positive A. The method of our proof does 
not allow to go further. Indeed, we have the following estimates, 

C{R) 



vr 

-+Ue- 



Ae 



R 



Ae 



ed^v" ( - + i*e - 



Ae 



L2 



A 



Therefore, one cannot hope that (|3.25p holds beyond t = ^ + i*e — Ae (with Ci 
proportional to A~^/^). On the other hand, if our final time is t — ^ + t^e — Ae 
with A > 0, we can prove the analogue of Prop. |3. 151 bv a "so long" argument (that 
is, p.25p with Ci proportional to A~^/^). As a result, we have the first part of 
Prop. 11.51 The last part follows from the remark we made above, that we know 5^ 



explicitly, therefore in particular its value at time i = + t*e — Ae. 



5. Anisotropic harmonic potential 
Consider the general harmonic potential in R", 

(5.1) V{x) = \{ujIx\ + ujIx\ + . . . + ujIxI), 

with LOj > for all j. It is isotropic when all the cjj's are equal, anisotropic 
otherwise. We suppose that the oj^'s take exactly d distinct values {2 < d < n), 
and renaming the space variables if necessary, we can assume that 

< LJl < UJ2 < ■ ■ ■ < ^d- 

We denote ij the multiplicity of Wj, 1 < j < d (ii + . . . + — n). At least two 
possibilities occur, as for the result one can hope for, corresponding either to Th. ll.2l 
or to Cor. 11.31 In the former case, one would be interested in the Cauchy problem 

, , f iedtu" + le^Au" = V{x)u' + e''" \u' l^" , {t, x)eR+x M", 

(5.2) <^ 2 

and in the latter, in 

ledtu" + -e^Au" = V{x)u' + e"'"|M^p'^M% (t, x) e 

(5.3) ' 2 



'"|t=o 



;»/2-^(e) ^ e"/2''' (e) ' 



where f,r'^^Il and — >0. We briefly discuss Eq. (|5.2p . and explain more 

precisely what happens for Eq. (|5.3p . 

For (I5.2|l . the same method as in Sect. [2]leads to the following phase, profile and 
operators. 



ip{t,x) 



1 " 

^9 Xl^J^i tan(wji), 



2 



(5.4) voit,x) = Yi 



Xi 



jJl yjcosiujjt) j \cos{uJit) ' ■ ' ■ ' cos(w„t) 

J'j{t) = sa\{ujjt) — icos{ujjt)dj, 
Hj{t) = uijXj cos{ujjt) + iesin{ujjt)dj 
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The first focusing occurs for t — ^5^; the solution focuses on the id-dimensional 
vector space defined by 

Ed = {xj = 0, Vj such that LOj = ujd\. 

Therefore, the critical index for the nonlinear term to be relevant in (|5.2p would 
be fc = dim = id- If > then the nonlinear term remains negligible up to 
time t = -TT— and before the next focusing, where the same discussion is valid. If 
k > max ij, then the nonlinear term is everywhere negligible, provided that no 

simultaneous focusings occur; indeed, the ujj part of the harmonic potential will 
cause focusing at times 

— H , K e Z. 

ZLUj UJj 

Two (or more) distinct Wj's can cause cumulated focusing if they are rationally 
related. To simplify the discussion, we now assume n = 2 and that uji and lu2 are 
irrationally related. In that case, focuses at time t = 1^ on the line {x2 = 0}. 
If fc = 1, then the nonlinear term becomes relevant near {{t,X2) = (2SjiO)}- The 
case of a focusing on a line was treated in [4] without potential, with an initial os- 
cillation that forces such a geometry for the caustic. With an anisotropic oscillator, 
the situation is technically much harder to handle. In |4j, no oscillation was present 
in the other space variable, and this variable could be considered as a parameter. 
In the present case, oscillations are always present in both space variables, so it 
is harder to measure the dependence of with respect to xi when it focuses on 
{x2 = 0}. We leave out the discussion at this stage. 

On the other hand, it is possible to understand (and prove) what happens for 
Eq. (|5.3p . Because we altered the time origin, the operators we now use write, 

Jj{t) = -j^ cos{ujjt) + i sm{ujjt)dj, 
Hj{t) = ujjXj sinlujjt) — ie cos{ujjt)dj. 

We also have an explicit formula for the linear solution (the analogue of Eq. (|1.3p ). 
which yields in particular Strichartz estimates. The solution of 



is given by 

n / \ 1/2 „ 

y^(t,x) = Y\ (—^ / e'^(*-^'^)/^/(2/)dy, 



where 



S{t, x,y)^}_^ ^7^^ [ costjjt - XjVj 

j=i J 



,2 



It is not hard to see that one can mimic the proof of Th. 11.21 to get the following. 

Theorem 5.1. Assume 2<n<5, ■i<(T< and let 2 < r < If n^2, 

there exists 6 > such that in either of the two cases, 

• a > cfq{2), or 

• \\fU<5, 
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the following holds (if 3 < n < 5, no additional assumption is needed). Denote 
■tpj- = Wj^^f, and 



2 



2 ^—^ taii{LOjt) 



Let be the solution of i5.3]) . Then for \t\ < ^ (that is, before refocusing) , and 



n , ^ 1/2 

- TT ( TT^ ) ^ 



e^o 11 V 2i7rsinaj,t / Vsinwii sinw„t 

j=i ^ J ' 



At time t — the solution focuses on iJ^, and the nonhnear term is neghgible 
(use the operators for all indexes j such that ujj = w^, and Jf, for the others). 
The nonlinear term will be relevant again only if there exists a time where the 
focusings caused by the different Wj's (1 < j < d) occur simultaneously, that is if 
there are positive integers ki, . . . ,Kd such that 

UJl "' LUd ' 

This means that the Wj's are pairwise rationally related. Therefore, at time t = ti, 
the caustic crossing will be described again by the scattering operator. Notice 
that when t approaches ti , the asymptotics given in Th. 15.11 has been modified in 
terms of Maslov indexes (for instance, since the crossing of Ed is linear, only linear 
phenomenon occur at leading order, that is precisely a phase shift measured by the 
Maslov index). More precisely, for (k^ — 1)11 /ud < t < ti, every ojj (1 < j < d) 
part of the harmonic potential has caused Kj — 1 (linear) caustic crossings, and 

uHt,x) - 1 I ^ , . r 



2i7r| sin[j,i| / Vsincjit sina;„t 



For ti <t < {Kd + l)TT/u!d, one has 



^ \2j7r| sinwjil y ysincjit' ' sina;„i 



^n.Xn 



J = 

and so on. 

If the Wj's are not pairwise rationally related, then only linear phenomena occur 
near caustics, and they are measured by Maslov indexes. 
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